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Phase transitions are usually treated as equilibrium phenomena, which yields telltale universality classes 
with scaling behavior of relaxation time and healing length. However, in second-order phase transitions 
relaxation time diverges near the critical point ("critical slowing down"). Therefore, every such transition 
traversed at a finite rate is a non-equilibrium process. Kibble-Zurek mechanism (KZM) captures this basic 
physics, predicting sizes of domains - fragments of broken symmetry - and the density of topological 
defects, long-lived relics of symmetry breaking that can survive long after the transition. To test KZM we 
simulate Bose-Einstein condensation in a ring using stochastic Gross-Pitaevskii equation and show that 
BEC formation can spontaneously generate quantized circulation of the newborn condensate. The 
magnitude of the resulting winding numbers and the time-lag of BEC density growth - both experimentally 
measurable - follow scalings predicted by KZM. Our results may also facilitate measuring the dynamical 
critical exponent for the BEC transition. 

Second order phase transitions are always associated with divergence of both the relaxation time ("critical 
slowing down") and healing length ("critical opalescence") at the critical point. These divergences imply 
inevitable non-adiabaticity when a system is driven across the transition and result in a finite size of 
domains that can coordinate symmetry breaking. That size is set by the healing length at the time when the 
system can no longer keep up with the externally imposed change of parameters that drive it through the 
transition, and adiabatic following gives way to impulse "freezeout" of its state. 

This mosaic of domains with independent choices of broken symmetry may result in topological defects. Their 
densities at formation will bear universal signatures of the underlying phase transition. Kibble-Zurek mechanism 
(KZM) is the theoretical framework that describes the dynamics of symmetry breaking in second order phase 
transitions by taking into account finite speed of propagation of the relevant information 1 " 5 . KZM leads to a 
quantitative estimate for the density of defects 3 " 6 . In particular, it predicts that their density should scale with the 
quench rate. The scaling exponent predicted by KZM is a function of the critical exponents of the underlying 
equilibrium phase transition. The basic ingredients for this prediction are the critical exponents - i.e., the 
universality class of the transition. Therefore, applicability of KZM spans physical phenomena of enormous variety 
and scales, starting from microscopic phenomenon of vortex formation in superfluid, right up to phenomena of 
astronomical scales, like formation of structures (cosmic strings, monopoles, etc.) in the early Universe. 

Kibble-Zurek mechanism has been tested in a number of primarily numerical studies, see, e.g., Refs. 7 " 14 . 
However, on the experimental side, while, as of now, data confirm key qualitative predictions of KZM (creation 
of topological excitations) 15 " 25 , its key quantitative prediction (scaling of their density with the rate of quench) has 
not yet been convincingly demonstrated (see, however, Refs. 23 " 25 , for suggestive indirect evidence). The difficulty 
involves controlling sufficient range of quench time scales as well as counting defects. 

In this article we pave a way around this long-standing hindrance. Within the setup of Bose-Einstein con- 
densation transition in an effectively one- dimensional ring, we demonstrate that KZM scaling laws can be 
observed involving stable and experimentally accessible quantities, namely, the variance of the winding number 
W and the non-adiabatic response time for density growth. Looking for signature of KZM in growth of bulk 
density is a novel perspective which allows bypassing the traditional difficulty of counting and resolving defects. 
Our demonstration is also a long overdue numerical verification of the original prediction proposed within the 
setup of superfluid circulation in an annulus 3 . 

Manifestations of symmetry breaking associated with BEC formation are spectacular and diverse 610 ' 12 ' 13 ' 1517 - 19 ' 26 . 
Formation of BEC breaks the (7(1) symmetry of the phase of the condensate wave-function. We show that when 
BEC forms in a ring cooled through the critical temperature, this symmetry breaking may result in spontaneous 
rotation in the newborn condensate. Persistence of the resulting flow is assured by topological stability of the 
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quantized winding number W. We show that Ws stabilize soon after 
the transition, and the variance of their distribution follows the scal- 
ing predicted by the KZM. Density growth of BEC also exhibits 
behavior consistent with KZM. 

Recently, persistent circulation of BEC in toroidal trap has been 
achieved experimentally by stirring the cloud 27 . In particular, with 
the advent of circular trapping potential for BEC 28 " 31 , the possible 
experimental testing of KZM is around the corner. Density growth in 
BEC formation with variable cooling rate has also been observed 15,32 . 
We show that a similar setup would allow both the testing of density 
growth scaling and determining the critical properties of the Bose- 
Einstein condensate 3 " 5 . 

Results 

The Model. We consider a BEC in a quasi- ID ring of circumference 
C, an idealization of quasi- ID toroidal geometry, e.g., se e 28 " 31,33 ' 34 . We 
model it using the stochastic Gross -Pitaevskii equation (SGPE) 35 " 37 : 



e( t )=-t/T Q 



(2) 



(i) 



where <j> = \(j){x)\e i0{x) is the condensate wavefunction and r\(x, t) is 
the thermal noise satisfying the fluctuation- dissipation relation (rj(x, 
t)rj*(x', t')) = 2yTS(x - x')3(t - t'), with y representing the 
dissipation, T the noise amplitude, g the non-linearity parameter 
and — e the chemical potential. We use dimensionless units given 
by H = 1, m= 1 (mass of a 87 Rb atom) and unit of time = l/co 0 with 
co 0 = 200 X 2ti Hz. We set T 0 =y _1 +y and £ 0 = 1 /v^ 7 ' 8 - 

Leaving aside the noise and dissipation, the above system can be 

described by the energy functional 8 = § . -\d x (j)\ 2 + C7(|0|) dx, 

I I 2 1 

where ^(|^|)= 6 |^| 2 + -g|^| 4 - Extremizing the energy functional 

we obtain <j> = 0 for e>0 and <j> = ^/|e|/gexp(/0) for e<0, where 6 
is the wave function phase (e = 0 is the critical point). We induce the 
transition by quenching e 



from an initial e > 0 to a final e < 0, and allow the system enough time 
to stabilize. The critical point is crossed at t = 0. Simulation of 
cooling of BEC by quenching chemical potential within the frame- 
work of SGPE has been shown to reproduce experimental results on 
defect (vortex) generation; see, e.g. Refs. 15 ' 35 ~ 37 . 

KZM in 1-D BEC Ring. When BEC is formed via cooling through 
the critical point, non-adiabaticity is enforced by diverging 
relaxation time t and healing length £ 3 " 5 : 

t = to/H vz ; { = {„/M v - (3) 

Here v and z are the critical exponents, £ 0 and t 0 are determined by 
the microscopic details of the system. According to KZM, as e 
approaches 0, after the instant —t the relaxation time t( — t) 
exceeds the timescale \e/^ t= _~ t of the change imposed by quench, 
and the state of the system freezes. Its order parameter behaves 
impulsively (i.e., remains effectively frozen) within an interval 
between + i and starts dynamical evolution again thereafter. KZM 
gives 

T (e(-f)) = |e/e| t= _-, = f- (4) 



For the linear quench in Eq. (2), we get from Eqs. (3), (4) 

1 = Z 0 \t Q /t 0 \^ 1 + VZ \ 



% U vz|l/(l + vz) 



(5) 



where \ = £ (i) . As e(t) becomes negative, condensate starts forming 
with a phase profile 0(x) consisting of random patches: phase is 
approximately uniform over the length-scale \. Within each such 
patch, its value is chosen independently and randomly (different 
stages of condensate formation in a uniform ring are shown in 
Fig. 1 upper row). Therefore, we can estimate the variance of the 
total phase 0q = § ring d0(x) within a torus of circumference C by 



W=1 





0 



-2 



Figure 1 | Condensate Growth in the ring. Sequences of isodensity surfaces for growing condensate for a single quench realization with t q = 0.01 for flat 
ring (top line, time flowing from (a) to (c)) and tilted ring with tilting angle a = 3° (bottom line, time flowing from from (d) to (f)). Color represents the 
phase of the condensate along the ring, showing magnitude of Wto be 1 and 2 respectively in the flat (c) and the tilted (f) rings. The 3D equivalent is 
constructed for a ring with isotropic harmonic transverse confinement with frequency co 0 = 200 X 2% Hz and other parameters as mentioned in the text 
for Figs. 2 and 3 (see "Methods" section for discussion on this equivalence). Formation of random equi-phased patches at the onset of condensation is 
clearly visible. 
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number. So, Eq. (5) predicts dispersion: 
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Figure 2 | Quenching into a BEC and its consequences. Time evolution of 
the winding number W, total number of BEC atoms N, total angular 
momentum L and specific angular momentum LI N for random instances of 
quench realizations with different t q . The quench is induced by ramping e 
linearly from -500 to +500 in each case (Eq. 2). For each t q , different 
quantities for the same realization are plotted in same color. The black 
vertical lines indicates t c = 0 (continuous, middle one) and the two on either 
sides (dashed) indicate + t y where t = sJtqTq, and t 0 = 7 -1 + 7- The second 
row shows the adiabatic- impulse transition in BEC growth. While t gives the 
KZM estimate for the impulse to adiabatic transition boundary, the instant 
denoted by A in the figure marks the point where N shoots up sharply. The 
thin blue lines indicate the equilibrium values of N for the corresponding 
instantaneous values of e. Prior to the quench, we thermalize the system at 
the initial value of the chemical potential. T = 10" 3 , y = 10" 2 and g = 0.05 are 
held constant for all time. 



considering the sum of ~ C j £ uniform random variables, each 

having the same variance n 2 /3 (0 taking any value between ±%). 
This implies Gaussian distribution for 0 C (in agreement with inset 

of Fig. 3a) with variance o 2 (0 c ) = (6 2 c ) = (n 2 /3) {cjt) when 

averaged over many realizations. As the wave function is single- 
valued, we must have 0 C = 2nW, where W is the integer winding 



a(W) = y/{m): 



1/8,-1/2 



L 0 



r 

so 



2x/3 



(6) 



using the mean field exponents z = 2 and v = 1/2 7 ' 8 . 

Spatial gradient of 0(x) gives local flow velocity. Therefore, W 
quantifies the net quantized circulation of the condensate around 
the ring. Thus, breaking of (7(1) symmetry leads to independent 
phase selection and results in a net superfluid circulation in the 
ring 3 5 . The scaling predicted above is valid when expected | W\ > 
1. It steepens for (\W\) ~ 1 or smaller; se e 23 ~ 25 ' 38,39 and references 
therein. 

We also observe that the scaling of i in Eq. (5) can be extracted 
from the growth of the condensate density. It shows a sharp "knee" 
behavior, revealing a characteristic response time as the transition 
point is crossed. The response time gives an effective i , which follows 
the scaling in Eq. (5). 

Simulation Results for Winding Numbers. We integrate Eq. (1) 
numerically to study the quench dynamics. Evolution of W, BEC 
particle number N, total angular momentum L and specific 
angular momentum L/N are illustrated in Fig. 2. We define total 
angular momentum L through single particle angular momentum 
operator. Assuming our wavefunction <j>(x) = \fN(j) single (x), me tota l 
angular momentunWonly the z-component is non-zero) is given 

by L = i£ t § c (l)*(x) — dx. Each quench in this figure continues 
from e = 500 to e = — 500, where it is held fixed for the rest of the 
time. Though the flow involves more mass as e decreases, W acquires 
a stable value right after the symmetry-breaking. Both L and N grow 
as long as e decreases, but L/N stabilizes to a steady value along with 
W. This is analogous to the build up of persistent current in 
superfluid He 3 " 5,40 . For spatially uniform density the simple 
relation L/N = W holds and we observe that at the tail-end of the 
quench (last row), when the density fluctuations are ironed out due 
to dissipation. Final kinetic energy of the BEC depends on the final 
steady value of W and scales with it quadratically. W stabilizes (first 
row) in the wake of the condensation, when N (second row) is still 
negligible. Thus, W retains phase information from soon after BEC 
formation at the time of symmetry breaking. Once stabilized, W is 
resilient to dissipation and typical ambient thermal fluctuations, due 
to its quantized nature. The scaling in Eq. (6) is verified by averaging 
\W\ over > 10 3 realizations. Here the quenches are also between 
e = + 500, except for the fastest one, for which it is done between 
e = ± 1000, so that the quench begins and ends outside the impulse 
region between + e. The simulation results compare favorably with 
the KZM prediction, as summarized in Fig. 3(a). A direct comparison 
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Figure 3 | Scaling Laws, (a) The scaling for (| W\) = ^/2/na(W) (averaged over > 10 3 realizations) with t q for C = 30 (open circle) and 120 (filled 
diamonds). The error bars indicate the standard deviations. Best fits (red line) in the Fig. give exponents 0.1262 ± 0.0026 for the C = 30 and 0.126 ± 
0.0075 for C= 120, compared to the expected v/2(l + vz) = 1/8 = 0.125 (Eq. 6). Random walk in phase prediction (| W\) - y/C, i.e., log«l Wl} c=120 ) - 
log((l W1)c=3o) ~ l°g2 (Eq.2), is also verified (the difference shown in the Fig. is ~ 0.308 compared to log2 ~ 0.301). The inset shows histograms for 
distribution of W. (b) ATvs t/t. The "knees", where AT shoots up (point A), show perfect overlap for different t q 's (values marked in the Fig.), and can be 
used for estimating critical exponents. All static parameters are same as in Fig. 2. 
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between the numerical results and KZM formula for our model (Eq. 
6) yields, o{ W)(KZM) = 5.94, 4.86 and 3.65 versus o{ W) (numerical) 
= 2.33, 1.83 and 1.35 for t q = 0.01, 0.05 and 0.5 respectively. With 
"naive" KZM (choices of phases over £-sized regions are completely 
independent), mismatch of this order is consistent with (actually, 
significantly less than) past observations 12,13 ' 41 . We shall comment 
on this further in the Discussion. The square-root dependence of 
(\W\) on C, Eq.(6) is also verified there, as the values of (\W\) 
obtained for C = 120 are approximately double of those obtained 
for C = 30 for the same t q s (Fig. 3a). 

Thorough initial thermalization is crucial in producing equili- 
brated initial state (and, hence, physically motivated final state). In 
particular, noise should randomize the phase along the entire cir- 
cumference C. Otherwise, winding numbers saturate (as we have 
observed with runs that did not attain this initial long-range therma- 
lization of the phase). This randomization of the phase is in effect a 
diffusive process that requires time ~C 2 , which renders accurate 
reproduction of the scaling behavior for large C numerically costly. 

Adiabatic-impulse transition is at the heart of the scaling laws 
predicted by KZM. Here we observe for the first time its manifesta- 
tion in BEC formation as a scaling-law for the length of the impulse 
reaction time (effective t; see, e.g. Ref. 41 ) for the condensate density 
growth (see however, similar plot for time-dependent spin alignment 
in spinor BEC 42 ). This is an important result, since density evolution 
in a condensation process is easier to study experimentally than the 
phase ordering associated with defect dynamics. We compare (Fig. 2, 
second row) the growth of N, with the instantaneous equilibrium 
value of N (thin blue straight-line segments), obtained from the 
relation e = —gN/C (valid for small g). Initially, e>0, and N is 
negligible. After crossing t c the instantaneous equilibrium value of 
N increases linearly with the rate N = C/giQ = 0.6 x 10 5 till t = 
500t q , where the ramp is stopped and equilibrium value for AT settles 
to 3 X 10 5 . But for the actual evolution, the length of the period from t 
= t c = 0 up to the point denoted by A (a sharp knee) in respective 
figures (see also Fig. 3b), should be proportional (and of the order of) 
i and - in our SGPE case - should scale as Tq 2 . Our results confirm 
this scaling, as shown by the overlap of N around the knee, when 
plotted against t/t for different Tq (Fig. 3b). After this knee point, N 
catches up rapidly with its equilibrium value. W and L/N also stabil- 
ize around A. 

Unlike other BEC relics of symmetry breaking (e.g., solitons), 
which are difficult to resolve due to thermal noise, and decay due 
to dissipation, W bears a very stable and readable signature of the 
underlying phase transition due to its topological stability and inte- 
ger nature. Statistics of W can be presumably studied, e.g. within the 
experimental setups such as 27 - 31 * 33 - 34 . Experimental study of the 
growth of BEC density with adjustable cooling rate parameter has 
been reported already by Esslinger group 32 (cooling by sudden ramp 
of confining RF field), and within toroidal geometry by Anderson 
group 15 (both with sudden and linear ramp of RF field). They 
observed the linear growth regime, as well as the "knee" feature that 
provides an experimental counterpart of the effective i. This should 
allow direct verification of the KZM scaling law and quantitative 
determination of the exponent for i. Moreover, experimental deter- 
mination of the dynamical exponent z may be possible employing the 
KZM formulae v/2(l + vz) and 1/(1 + vz) for the measured values of 
(| W\) and effective t exponents, Fig. 3(a) and 3(b), respectively, and 
solving for v, z. The scaling law involving AT may be amenable to more 
accurate experimental determination, since the exponent for t scal- 
ing (i ~ T Q 2 ^j is larger than that for the scaling of (a(W)~ Tq . 
Note that for real BEC transition in 3D, theoretical prediction is v = 
2/3 (close to experimental value, 0.67 ± 0.13, 44 ) and z = 3/2. Finally, 
the sample to sample fluctuations in the growth of N are much 
smaller than those of W (see Fig. 2), and hence requires less averaging 
over realizations. 



Discussion 

Our simulations reproduce scalings of the winding number and i 
predicted by KZM, but the "naive" random walk argument including 
Eq. (6) we have used to estimate < | W\ > results in an overestimate 
by a factor / ~ 2.5. This discrepancy deserves a comment both 
because it is there and because it is surprisingly small: Factors / ~ 
10 were reported in the past 12,13,41 . 

We first note that KZM estimate - one f - sized defect fragment per 
domain of size £ - is general, and this generality is obtained at the 
price of focusing on what is dominant (scalings) and ignoring (sub- 
dominant) microphysics. For instance, it is clear that the estimate of 
£ will depend on the noise temperature and on the nonlinearity 
parameter. Moreover, density of defects will depend on their nature, 
and there are systems (such as 3 He) where there are many different 
types of defects that can be created by the same transition, and there 
is no reason for their densities to be the same. 

Detailed studies show that such additional inputs affect density 
estimates but are either subdominant (e.g., there is a logarithmic 
dependence on noise and nonlinearity, as can be expected from the 
nature of the instability after the quench 43 ) or are too complicated for 
precise treatment ( 3 He). However, these additional inputs do not 
change scaling of f , and therefore, scaling of the number of topo- 
logical defects with the quench rate. 

This focus on scalings is very much in the spirit of the renorma- 
lization approach to second order phase transitions: Microphysics 
sets dimensionfull inputs that determine healing length and relaxa- 
tion time, but scalings of these two quantities near the critical point 
are independent of microphysical details. This allows for the 
classification of second order phase transitions in terms of critical 
exponents. Similarly, KZM predicts scalings, but only estimates 
prefactors. 

In contrast to past numerical experiments (where defect separa- 
tions were typically/^, with/— 10) we have seen KZM estimates of 
winding numbers only a factor of about 2.5 larger than the numerical 
estimates. This improved accuracy is because - in contrast to past 
exercises of this sort - in the constrained problem we have addressed 
we were able to be more precise in taking "SGPE microphysics" into 
account. Thus, we have estimated that domains of size £ will have 
phases that differ by n/ ^/3, and the geometry of the problem allowed 
us to avoid some of the difficult questions concerning dynamics of 
the post- quench phase ordering. However, these inputs are still only 
estimates: It is clear from Fig. 1 that different "beads" of the newly 
formed BEC do not have the same size. Moreover, it is quite likely 
that early on - around the time when the condensate begins to 
encompass the whole torus - winding number can change (see 
Fig. 2). Last not least, the instant when this seems to happen scales 
in accord - but does not coincide - with the estimate of i. 




0.5 1.0 

a(7t/180) 

Figure 4 | Winding up a tilted BEC ring, (a) ( I W1 ) vs a for t q = 0. 1 (other 
parameters same as in Fig. 2). (b) Suppression of scaling behavior of ( I W\ ) 
with Tq for a = 1°. The inset shows scaling behavior for smaller tilt of a = 
0.5° with a larger exponent 0.158 ± 0.011. 
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We conclude this part of the discussion by noting that KZM is a 
general theory of the dynamics of second order phase transitions, but 
that - like the universality classes on which it is based - it largely 
ignores microphysical details and does not aim to predict dimen- 
sionfull prefactors. Thus, while a more careful estimate of the 
expected magnitudes of winding numbers are possible in the simple 
model we have examined, the real aim of KZM and a good way to test 
it is to verify that it correctly predicts scalings of the topological relics 
of the transition from the critical exponents (which, for the real BEC, 
are anyway different than in the case of our SGPE model). Toroidal 
BEC traps are needed for this, and they are not the "popular model". 
Fortunately, our study suggests an alternative, more limited but also 
potentially more accessible possibility - testing for the scaling of 
delay in BEC density growth. 

Predictions of KZM can be substantially perturbed by natural 
sources of experimental imperfections in the form of inhomogeneity 
in the potential. An overall non-uniformity may have many causes. 
To be specific, we shall think of tilting the torus in the gravitational 
field. This is represented by an additional potential of the form 
V(S) = mgsm(oc)R[l — cos(«9)] where a represents the tilting angle, 
8 = x/R is the angular coordinate denoting the position on the ring, R 
= Clin and g the gravitational acceleration. Different stages of con- 
densate formation in the tilted ring are shown in Fig. 1 (lower row). 
Fig. 4(a) shows suppression of W as the function of a. Inset of 
Fig. 4(b) shows that the scaling behavior still persists for very small 
tilting of a = 0.5°, but the exponent increases. The main Fig. 4(b) 
shows the disappearance of the KZM scaling for slightly bigger tilting 
(a = 1°). Suppression of W could be caused in the tilted ring due to 
the competition between the finite velocity v F of the critical front 
(determined by the critical condition set locally by the inhomogen- 
eity 6 ) and the relevant sound velocity v = \ / t (t) at which the cor- 
relation is established between the condensate phase. Formation of 
phase-inhomogeneity (local flow) is suppressed wherever v>v F , 
since phase correlation is maintained throughout the condensation 
process there: the newborn condensate selects the same phase as the 
existing condensate due to local free- energy minimization 6,45 ' 46 . This 
implies steeper fall of (\W\) with t q 6 . But in our case, the above 
condition is not met for a = 0.5°. The same criteria also cannot 
explain the suppression of W for a = 1° (Fig. 4b). One possible 
explanation of this suppression may lie in the dissipation of kinetic 
energy prior to the formation of the complete BEC ring. The local 
flows generated by KZM at the bottom of the ring are susceptible to 
dissipation before the topological constraint (that protects the quan- 
tized circulation) is imposed, i.e., before the BEC is formed comple- 
tely up to the ring top. Such early energy loss might leave the 
condensate with kinetic energy insufficient for quantized circulation. 
We note that tilting a BEC ring experimentally and applying that to 
offset inadvertent horizontal misalignment has been achieved 
recently 34 , so while we have not developed the theory of various 
imperfections (exemplified here by tilting in the gravitational field), 
control with the accuracies that may compensate for (or avoid) such 
problems seems possible. 

For very slow quenches (large t q ), however, | W\ doesn't vanish in 
the tilted case, as the density doesn't grow on the top till the critical 
front reaches the top, and topological protection doesn't apply till then. 
If this period is long enough, thermal fluctuations may create sufficient 
random walks in phase to induce a net circulation by the time the 
density becomes significant everywhere. The resulting circulation may 
hence even grow with Tq in this regime, as seen from the Fig. 4. 

To summarize, we showed that temperature quench through the 
critical point can produce spontaneous circulation of BEC in a ring 
and verified long standing scaling predictions of KZM 3 . Our demon- 
stration involving winding number and condensate density paves a 
way for experimental verification of KZM scaling laws. It also pro- 
vides prescription for experimental determination of the critical 
exponents of the underlying BEC transition. 



Methods 

We simulate the quasi- ID BEC by integrating 1-D Stochastic Gross Pitaevskii 
equation (Eq. 1) numerically using 4th order adaptive step-size Runge-Kutta algo- 
rithm for > 1000 noise realisations. The observables are extracted through averaging 
over the ensemble 15 ' 50 . The initial state for the evolution of Eq. (1) is a thermal state 
obtained integrating Eq. (1) with constant e fixed to its large initial value, for an 
integration time > 1000t 0 to ensure thermalisation. The temperature quench of the 
thermal cloud is simulated through a ramp in the chemical potential ( — e), keeping 
the noise amplitude T and the damping y constant. This approximation captures the 
essential universal aspects of the dynamics. The symmetry breaking underlying the 
transition is induced by the change in e (see the discussion following Eq. 1). Changes 
in T and y give only non-universal quantitative corrections 713 . A similar method was 
successfully used in Ref. 15 in the framework of the stochastic (projected) Gross- 
Pitaevskii equation (SPGPE) to model evaporative cooling of a pancake shaped 
thermal cloud and study the statistics of trapped vortices and found good agreement 
between experimental and numerical results. The SGPE and the SPGPE have been 
used before to model experimental results with success 1535 ' 36,47 " 49 . 

To make direct connection of our simulation with experiments, we note that the 
one dimensional approximation of a thin 3D torus with thickness d(d<^£) is obtained 
by integrating out the transverse (y and z in our case) degrees of freedom. Such a 
confinement- induced quasi-one dimensional Bose system is equivalent to a purely 
one-dimensional system with a modified interaction strength g given by the con- 
finement frequencies. For isotropic harmonic confinement in the transverse direction 
with confinement frequency co 0 we have g = 2% 2 aj m a\, where m is the mass of the 
atom, a is the scattering length, and a§ = ^ti/mco 0 . The three dimensional wave 
function (f>(x, y, z) may be reconstructed from the one-dimensional wave function 
(j)(x) by assuming Gaussian profile in the transverse directions: (j)(x,y, z) = 0(x)0o(y» 
z), where (j) 0 (y,z) = e^ 2 +z2 )/ 2a l . Such reconstructions are shown in Fig. 1 and in 
the Supplementary Movie. The Supplementary Movie shows the real-time dynamics 
of the condensation process illustrating the key features of the dynamics germane to 
verification of KZM and determination of critical exponents, namely, stabilization of 
W at the wake of the condensation and the knee-like behavior of the density growth 
(also clearly shown in the corresponding Supplementary Fig.). The simulation is done 
by integrating dimensionful SGPE with parameter values: 

a 0 = 100 X a B (a B = Bohr radius), m = mass of a Rb atom, co 0 = 200 X 2n, y = 0.01 
and T = 10 nK. All these parameters are compatible with the experimental set up in 28,29 . 
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